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Abstract. We summarize the existing experimental data on electromagnons in multiferroic 
/?Mn205 compounds, where R denotes a rare earth ion, Y or Bi, and discuss a realistic 
microscopic model of these materials based on assumption that the microscopic mechanism 
of magnetically-induced ferroelectricity and electromagnon absorption relies entirely on the 
isotropic Heisenberg exchange and magnetostrictive coupling of spins to a polar lattice mode 
and does not involve relativistic effects. This model explains many magnetic and optical 
properties of /?Mn205 manganites, such as the spin re-orientation transition, magnetically- 
induced polarisation, appearance of the electromagnon peak in the non-collinear spin state 
and the polarisation of light for which this peak is observed. We compare experimental and 
theoretical results on electromagnons in /?Mn20s and 7?Mn03 compounds. 
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1. Introduction 

Multiferroic materials that exhibit simultaneous magnetic and ferroelectric order have 
attracted much attention recently because of the fundamental interest of systems with coupled 
order parameters and because of their potential for cross electric and magnetic functionality 
EH SIS migl 12. Many recently discovered multiferroic s, e.g. TbMn0 3 QD, TbMn 2 5 [U, 
and N13V2O8 ifTOl . are improper ferroelectrics in which the electric polarisation is induced 
by spin ordering. These materials show striking cross-coupling effects such as magnetic field 
induced polarisation switching flUIH and giant magnetocapacitance ifTTl and a strong coupling 
between spin and lattice excitations that leads to electric dipole excitation of magnons, 
or electromagnons 01 [12l [131 EHL which is the subject of this paper. The fundamental 
interest in multiferrocity also derives from the strong interplay between magnetic frustration, 
ferroelectric order and unusual symmetry breaking in phase transformations that characterize 
these materials ff6l[T5i 

The known microscopic mechanisms of magnetically-induced ferroelectricity include 
lattice distortion (exchange striction) and redistribution of electron density in response to spin 
ordering. Such processes occur locally in all magnetic materials. However, only when a spin 
ordering breaks inversion symmetry do these local electric dipoles add into a macroscopic 
electric polarisation. Spin orders that break inversion symmetry are rare and the best systems 
to look for them are frustrated magnets, where competing interactions and the geometry 
of spin lattice favor unconventional magnetic states. In most of the recently discovered 
multiferroic materials, such as TbMn03, Ni3V20s, MnWCU and CuO, competing interactions 
force spins to form a cycloidal spiral. This non-collinear spin order breaks inversion symmetry 
and activates antisymmetric Dzyaloshinskii-Moriya interaction proportional to the cross- 
product of spins, S\ x S2 lfT6l . The concomitant lattice and electronic distortion induces 
electric polarisation lfT71[T8l[T9Tl . 

When inversion symmetry is broken by a collinear magnetic ordering the strongest spin 
interaction that can shift ions and polarise electronic clouds is the symmetric Heisenberg 
exchange, proportional to the scalar product of spins, S\ ■ S2. This mechanism was proposed 
to explain multiferroic properties of i?Mn205, where R denotes a rare earth ion, Y or Bi, and 
orthorombic manganites showing the E-type antiferromagnetic ordering ll2"0ll2Tll2"2~l| . 

These two microscopic mechanisms of magnetically-induced ferroelectricity give 
rise to two different forms of phenomenological magnetoelectric coupling: electric 
polarisation induced by a low-pitch spiral is described by the third-order coupling term 
P(L\dL2 —LidLi), where P is electric polarisation and L\ 2 are magnetic order parameters 
describing the sinusoidal and cosinusoidal components of the spiral ||23l I241 . while the 
coupling working in collinear spin states has the form P^L^—L^), where L\ and L2 are 
components of a two-dimensional irreducible representation describing magnetic states with 
opposite electric polarisations [|25l[20ll22l| . 

The magnetoelectric interactions that induce electric polarisation in magnets can also 
couple oscillations of magnetisation to polar lattice vibrations. The oscillations of polarisation 
at magnon frequency and vice versa give rise to dynamic magnetoelectric effects, such as 
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Figure 1. (Colour online) Feynman diagrams describing photo-excitation of magnons by the 
electric field of light (here photon and phonon are represented, respectively, by dashed and 
solid lines, while the wavy line corresponds to magnon): (a) Photoexcitation of two magnons 
and one phonon via the fourth-order magnetoelectric coupling (the Lorenzana-Sawatzky 
mechanism), (b) Excitation of the two-magnon continuum via the third-order magnetoelectric 
coupling ('charged magnons') and (c) Photoexcitation of a single magnon (electromagnon) 
due to the third-order magnetoelectric coupling, where photon with the frequency co and zero 
wave vector scatters off the static spin modulation with the wave vector Q producing a magnon 
with the same wave vector and the frequency (0. 



electromagnon excitations. In usual magnets, an oscillating electric field of photons can 
excite a three-particle continuum consisting of two magnons and one phonon 11261 . This 
process results from the fourth-order spin-lattice coupling (see figure [Tk). The third-order 
couplings in multiferroics, discussed above, make possible photo-excitation of two-magnon 
continuum without a phonon ('charged magnons' 112710 shown in figure [TJ)). Replacing one 
of the magnons by the static modulation of spin density appearing in the ordered spin state, 
we obtain a process that converts photon into a single magnon, which is the electromagnon 
(see figure [Tfc). This process is usually mediated by a polar phonon linearly coupled to both 
magnons and light, which for low-frequency phonons can lead to a resonant enhancement of 
the photo-excitation of electromagnons. 

As there are two possible contributions to the polarisation — from ionic displacements 
and from electronic density redistribution — one can think of two corresponding electric 
contributions to electromagnon — from phonons and from electronic excitations. This means 
a transfer of the electric dipole spectral weight from phonons (hco ~ 10 — 100 meV) and/or 
electronic excitations (hco ~ 2 eV) down to magnons (hco ~ 1 — 5 meV). Such a transfer, in 
turn, leads to a step-like anomaly in the temperature dependence of the dielectric constant. 
Current research on electromagnons, including this paper, is focused mainly on magnon- 
phonon coupling while magnon-electron aspect is much less explored as it is generally 
expected to be weak because of the large energy difference between the electronic excitations 
and the magnetic excitations. 

The possibility of electromagnon excitations in multiferroics has long been anticipated 
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theoretically [1], but only recently were they observed in experiment. Pimenov et al lfT2l 
reported observation of electric active modes at magnon frequencies in GdMn03 and TbMn03 
which exist only in some magnetically ordered phases and can be suppressed by magnetic 
field. Further, they found in GdMn03 the spectral weight transfer from the lowest frequency 
phonon down to the electromagnon mode ll28ll . Sushkov et al |[T3l reported observation 
of electromagnons in YM^Os and TbM^Os. Electromagnons in both these 7?Mn205 
compounds have very similar spectrum and exist only in magnetic ferroelectric phases that 
proves their magnetic origin. Valdes Aguilar et al ll29ll showed that electromagnon absorption 
in Euo.75Y .25Mn03 occurs over a broad band with two peaks, both of which exist only in the 
magnetic ferroelectric phase below 30 K. Pimenov et al |[30l explored a composition set 
Eui A YJVIn03 for < x < 0.5 and confirmed main features of electromagnons in 7?Mn03. 
Kida et al OTTl found out that the electromagnon polarisation in DyMn03 stays with the lattice 
when the spin plane is rotated by the external magnetic field. Electromagnons, observed so 
far only in non-collinear spin phases of 7?Mn03 and 7?Mn205 compounds, have common 
features for both families: they are active only in one polarisation e\\a axis for 7?Mn03 and 
e\\b axis for i?Mn205, where e is electric field of light; well defined peaks exist only in the 
low-temperature magnetic ferroelectric phase. A theory of electromagnons for the circular 
magnetic spiral was developed by Katsura et al lfl4l and its continuum version was given by 
de Sousa et al 0211331 who applied it to describe low-energy excitations in BiFe03. Fang and 
Hu [34] discussed electromagnon absorption in 7?Mn205 compounds assuming Heisenberg 
interactions between spins. 

The outstanding fundamental questions for electromagnon are the microscopic origin of 
these novel excitations (Heisenberg or Dzyaloshinskii-Moriya type exchange) in the different 
classes of compounds (i?Mn03 and 7?Mn205), the explanation of the observed selection rules, 
and whether electromagnons may occur in a wider range of materials. The practical issues are 
enhancing the magneto-capacitance effect and its temperature range and possibly applying 
these new excitations to metamaterials and/or achieving negative index of refraction in the 
magnon range of frequencies. 

The remainder of this paper is organized as follows. In section [2l we discuss optical 
absorption spectra of multiferroic YM^Os and TbM^Os compounds taken at different 
temperatures. We present experimental evidence allowing us to identify the low-frequency 
peaks appearing in the non-collinear phase of these materials as electromagnons. For 
comparison, in section [3] we present data on electromagnons peaks in the spiral multiferroic 
material Euo.75Yo.25Mn03. In sectionHl we discuss a simple microscopic model of YM^Os 
manganites, based on isotropic Heisenberg exchange interactions between spins, with which 
we describe magnetic, multiferroic and optical properties of these compounds. We also briefly 
discuss phenomenological description of magnetic orders and magnetoelectric coupling in 
7?Mn205 compounds based on their symmetry. In section[51 we discuss the microscopic origin 
of the electromagnon peaks in both i?Mn205 and 7?MnC>3 and how the electromagnon relates 
to the spontaneous polarisation. Finally, we summarize our experimental and theoretical 
results. 
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Figure 2. (Colour online) Optical conductivity of YMr^Os and TbMn?Os for the electric field 
of light e\\b in three phases. Strong peaks at 113 and 97 cm -1 are the lowest phonons, other 
peaks are electromagnons. 



2. Far infrared spectroscopy of i?Mn20s compounds 

The family of i^Mr^Os compounds has long been studied ll35l . Small variations of the 
exchange integrals with temperature lead to the complex phase diagram for this multiferroic 
family 11361 . Earlier spectroscopic works revealed far infrared absorption modes activated at 
low temperatures in EuMn 2 5 (371, YMn 2 5 flUD and GdMn 2 5 [391. Recently, we have 
shown that in YM^Os and TbM^Os strong electric dipole active modes emerge at magnon 
frequencies in the lowest in temperature ferroelectric phase lfP3~l . As in the work by Pimenov 
et al 021 on 7?Mn03 compounds, we assign this modes as electromagnons. 

In figure |2l we show the optical conductivity spectra of YM^Os and TbMmOs at 
three different temperatures corresponding to three magnetic/ferroelectric phases. These 
spectra were obtained by fitting measured transmission spectra [TT3l with the model dielectric 
function Il40l . assuming all modes to be of the electric dipole nature. The spectra measured at 
7 K (blue curves) show the characteristic strong and sharp (electromagnon) absorption peaks 
at lowest frequencies. In this phase, the spontaneous electric polarisation is relatively small 
and the angles between neighbouring spins along the b axis are large. The spectra measured 
at 25 K (green curves) have one broad absorption peak near 20 cm -1 . In this phase spins 
are almost collinear and electric polarisation is large. Red curves, taken just above the Neel 
temperature, show a single broad absorption band below the phonon frequencies. 

Identifying the low-frequency excitations as electromagnons requires addressing several 
questions. First, to avoid confusion with possible transitions between /-levels of rare 
earth ions, we have studied YM^Os. The second issue is electric versus magnetic dipole 
(antiferromagnetic resonance) activity. We measured transmission spectra for various mutual 
orientations of the electric and magnetic field of light (respectively, e and h) with respect to 
the crystal axes and we found the absorption only for e||6||P, where P is the spontaneous 
polarisation vector, independent of the orientation of h, which implies that the excitations are 
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Figure 3. (Colour online) Dielectric constant of TbMmOs from fits of infrared spectra 
(lower curve) in comparison with kHz measurements. The whole step-like anomaly is due 
to electromagnons (figure[2k, blue curve). 



electric dipole active. Can these resonances be new phonons, activated in the low temperature 
phase? We have performed shell model calculations that put the lowest phonon near 100 cm -1 
for any reasonable parameters. Thus, we believe that we can reliably identify these low- 
frequency peaks as electromagnons. 

Another check of the electromagnon origin of this modes is comparison of their 
contribution to the step-like anomaly in the temperature dependence of £\, calculated using 
the Kramers-Kronig relation, with the measured one. Figure [3] shows such a comparison. We 
have chosen TblVh^Os for this purpose because of the larger sample size and higher frequency 
of the electromagnon peaks. It is clear from figure[3]that the whole step-like anomaly in £\ (T) 
comes entirely from the sharp electromagnon peak. 

The frequency and temperature behaviour, presented in figures [2] and [3] is typical for a 
set of i?Mn 2 5 compounds: R = Er, Ho, Y, Dy, Tb, Gd and Eu El 11 [II [33 [33 11. ThC 
appearance of electromagnon peaks as well as the £ step-like anomaly seem to be correlated 
with the transition into the non-collinear spin state. Notably, BiM^Os whose spin ordering is 
nearly collinear at all temperatures, shows neither the £ step-like anomaly nor electromagnon 
absorption H3l . 

An inelastic neutron scattering study of YM^Os is in progress. Preliminary data by 
S.-H. Lee et al P4l show several scattering peaks in energy scans at the wave vector of the 
static spin structure. A strong neutron feature is observed at 1 meV in good agreement with 
the sharp low frequency feature in the low temperature infrared spectrum (figure [2k). 

Katsura et al lfl4l predict that the electromagnon originating from antisymmetric 
exchange has a 'transversal' polarisation with respect to the spontaneous polarisation: e _L 
P. The observed polarisation selection rule for electromagnons in 7?Mn205 compounds is 
'longitudinal': e\\P. We will show that such 'longitudinal' electromagnon can be obtained in 
a model based on the isotropic Heisenberg exchange. 
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Figure 4. (Colour online) Optical conductivity of Euo.75Yo.2sMn03 for e\\a in the 
ferroelectric (blue), the spin density wave (dashed), and the paramagnetic (red) phases. 
Electromagnon band consists of a broad background and two peaks. 



3. Far infrared spectroscopy of i?MnC>3 

It is interesting to compare electromagnons in the two families of multiferroic manganites. 
We begin with Euq .75 Yq 2sMn03 — a compound that mimics lattice parameters of TbMn03 
but does not have /-/ transitions in the far infrared [|29l . To extract the parameters of the 
oscillators, we fit the transmission spectra with a Lorentzian model of the dielectric constant 
e(co) for electric dipole transitions: 

e{a>) = Co, + £ — — -{ (1) 

where £«, is the high frequency dielectric constant, j enumerates the oscillators, Sj is the 
spectral weight, 0)j is the resonance frequency, and jj is the damping rate. 

Figure |4] shows optical conductivity spectra for three phases of Euo.75Yo.2sMn03. 
Figure |5] shows temperature dependence of the fit parameter Sj (spectral weight of the peaks 
in figure HI). A broad absorption (electromagnon) band exists well above the Neel temperature 
47 K (red curve in figure|4]and curve 3 in figure|5]). Upon cooling, the spectral weight S of this 
broad band is growing down to 7>£=30 K and stays at this level at lower temperatures. Such 
a behaviour of the background can be seen from the spectra — the absorption at 40 cm , 
the frequency least affected by two absorption peaks, stays constant at all T < Tfe- This 
absorption produces smooth growth of the e a (T) at Tfe < T. T^=41 K is the inflection point 
of the curve 3 which is an evidence for the magnetic origin of this low frequency broad electric 
absorption. Two electromagnon peaks emerge sharply at the Tfe (blue curve in figure |4] and 
curves 1 and 4 in figure [5]). They produce all further growth of the £ a (T) at T < Tfe- The 
electromagnon spectrum of TbMn03 is very similar to that one of Euo.75Yo.25Mn03, except 
the large electromagnon peak is at 60 instead of 80 cm -1 [|43l . 

It is also interesting to follow the temperature dependence of the redistribution of the 
electric dipole spectral weight S. Since electromagnons result from a small admixture of 
phonons to magnons, the total spectral weight should be conserved. Comparing the curves 2 



Electromagnons in multiferroic RMn20$ 



8 





i i „— 1.38 






' T T " 

FE IN ? 

■ 5 ; : f 

■ II! 2 




1x10 


■ I i O 1.36 

" i i I 


■ 


CM 

E 


V ! 6 50 100 150 200 250 
i Temperature --- ----^'"^ 

4 l i ^ — 

'\Ji Eu Y MnO 

■^w^^ak ! 0.75 0.25 3 


5x1 3 


-3 \ A\ 

■ y\ 


1 — ■— peak 25 cm" 1 

2 — ■— phonon 1 20 cm" 1 






3 —A— background 70 cm 1 

4 — *— peak 80 cm" 1 







5 — ■— 1+3 + 4 

— r- , A 



50 100 150 

Temperature (K) 



Figure 5. (Colour online) Spectral weight of the absorption peaks below 140 cm -1 in figure|4] 
The numbers from 1 to 5 enumerate the curves. Frequencies in the legend are the lowest 
temperature resonance frequencies of each mode. Inset: Total spectral weight of the eight 
phonons above 140 cm -1 . 

and 5 in figure [51 one can estimate how much spectral weight the 120 cm -1 phonon is loosing 
and how much electromagnons acquire. As electromagnons gain more spectral weight than 
the phonon is loosing, we checked the rest of the phonon peaks for the same polarisation of 
light. The inset in figure[5]shows that, indeed, phonons lose just enough of the spectral weight 
to conserve total spectral weight. 

In their recent inelastic neutron scattering work, Senff et al ll45l reported a set of modes 
at the incommensurate zone center of TbMn03. The frequency of the lowest mode is equal to 
the frequency of the lowest infrared peak (24 cm -1 ). We agree with Senff et al in assignment 
of this mode as electromagnon. However, despite the satisfied polarisation selection rule 
e _L P, this is not the electromagnon predicted by Katsura et al flU, as we discuss below. 

4. The model of 7?Mn 2 5 

The presence of two different types of magnetic ions and geometric frustration of 
spin interactions give rise to rather complex magnetic structures 7?Mn205 compounds. 
Nonetheless, a number of salient properties of these materials, such as the magnetically- 
induced electric polarisation, photo-excitation of magnons as well as the spin re-orientation 
transition, can be understood within a simplified microscopic model, which we discuss in 
this section. Our starting point is the assumption that multiferroic and optical properties of 
these materials are governed by the isotropic Heisenberg exchange, although we do include 
magnetic anisotropy to explain the spin re-orientation transition that has a strong effect on 
the low-frequency absorption spectrum. We discuss ordered spin states of the model, the 
mechanism of magnetoelectric coupling and calculate the optical absorption spectrum at 
magnon frequencies for different magnetic states. 
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We start by considering a single magnetic ab-layer including Mn 3+ and Mn 4+ ions. The 
model describes interactions between the spins and their coupling to a polar phonon mode: 

H = jY,i,jJij (P) {&i ' Sj) ~ ^HiaKia (Si ' ha) ~ 

-Ziiii{Si-H)+v(^-PE-^y (2) 

Here the first term is the isotropic Heisenberg spin exchange, while the second term is the 
single-ion anisotropy. The antisymmetric Dzyaloshinskii-Moriya exchange as well as other 
types of anisotropic exchange interactions are not included, as in the scenario discussed below 
they play no role. We assume for simplicity that in the ordered state all spins lie in the ab 
plane [1211 and neglect the small out-of-plane components found in recent neutron diffraction 
experiments on single crystals ll46l l47l . Thus, the easy and intermediate magnetic axes on 
each Mn site (a = 1,2) lie in the ab plane, while the hard axis kn \\ c. The third term in 
© is the interaction of spins with an applied magnetic field and the last term describes the 
dielectric response of the system, where is the 'bare' dielectric susceptibility related to 
the polar lattice mode (not including the magnetic contribution calculated below) and X2 is the 
remaining part of the dielectric susceptibility of non-magnetic origin. Finally, V is the system 
volume. 

The coupling between the spins and the polar phonon mode results from the dependence 
of the exchange coupling on the electric polarisation, which in i?Mn205 is parallel to the b 
axis: 

Jij(P b ) = JijiO) +J' ij (0)P h + \j'^)Pl + ... (3) 
The last two terms give rise to the cubic and the quartic magneto-electric couplings. 

4.1. Magnetic ordering and spin re-orientation transition 

Here, we adopt the model of Chapon et al lt2TI with 5 exchange constants between pairs 
of nearest-neighbour Mn ions: J\ and J 2 couple Mn 4+ ions along the c direction, Jj and J$ 
couple the spins of neighbouring Mn 3+ and Mn 4+ ions and J4 is the coupling between two 
neighbouring Mn 3+ ions (see figure©. 

Figure [6] shows the minimal energy spin configurations obtained by the numerical 
minimization of the spin energy © on the subspace of the commensurate magnetic states with 
the wave vector Q = (1/2,0,0) for two different sets of exchange constants. The exchange 
constants J4 and ^5 were chosen to be positive and large compared to other exchange constants, 
which gives rise to antiferromagnetic zig-zag chains along the a axis with nearly collinear 
spins (marked by dashed lines) observed in neutron experiments ll2Tll47l . 

The angle between spins in neighbouring chains sensitively depends on the ratio between 
the interchain coupling J3 and the magnetic anisotropy parameters Kj. We assume that the 
easy magnetic axis is parallel to the a axis. The interplay between magnetic anisotropy and 
interchain interaction determines the angle between spins in neighbouring antiferromagnetic 
chains. 
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Figure 6. (Colour online) Minimal energy spin configurations for J4 = J5 = 40 K and the 
anisotropy parameter I,,(Mn 3+ ) = 0.6 K for all Mn 3+ ions (red) and K a (Mn 4+ ) = 0.1 K for 
all Mn 4+ ions (blue). The value of the interchain coupling ^3 is —2 K for the structure in panel 
(a) and —4 K for the one panel (b). 

We first note that if spins in each antiferromagnetic chain would be perfectly collinear 
than the interchain interactions would cancel as a result of geometric frustration. Conversely, 
the interchain coupling 73 results in spin rotations, which destroy the collinearity of spins 
in each chain. Consider, for example, the spin of the Mn 4+ ion, marked in figure [6£a) by an 
arrow. Because of a nonzero angle between spins in the neighbouring a-chains, the interaction 
of the marked spin with the spin of the Mn 3+ ion from the neighbouring chain will give rise 
to a rotation of the marked spin. These small spin rotations lift the frustration and lead to 
some energy gain due to interchain interactions. It easy to show that this energy gain is 
maximal when spins in neighbouring antiferromagnetic chains are orthogonal to each other. 
Thus, while the magnetic anisotropy favours an almost collinear spin configuration, interchain 
interactions favour the 90° angle between spins of neighbouring chains. 

This competition gives rise to a very strong sensitivity of the angle between spins in 
neighbouring chains to the interchain coupling J3. For weak interchain coupling J3 = —2 K, 
this angle is relatively small and minimal-energy spin configuration shown in figure [6£a) is 
similar to the one observed in the high-temperature 'collinear' phase of YM^Os by Chapon 
et al [|2T|. A small change in 73 from — 2 K to —4 K transforms the configuration shown 
in figure [(Ja) into the one shown in figure [6tb), which may explain the spin re-orientation 
observed in i^M^Os with R =Tb, Ho, Dy and Y, provided that the interchain coupling is 
temperature-dependent. Although the rotations that make spins in each chain non-collinear 
are barely visible, they are sufficient to produce the large changes in the spin configuration, 
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since the magnetic anisotropy is relatively weak. In 7?Mn205, the spin re-orientation transition 
is accompanied by the loss of commensurability of the spin structure in the a and c directions 
(which is also a consequence of magnetic frustration). This latter aspect of the transition 
is, however, less important for the photo-excitation of magnons, discussed below, than the 
re-orientation of spins. 

4.2. Magnetically -induced polarisation 

Minimizing © with respect to P y , we obtain expression for the magnetically-induced electric 
polarisation, 

y(°) 

n*-%-JLW(StSj), (4) 
which only involves scalar products of spins. 




Figure 7. A cartoon of the magnetic ordering in the high-temperature collinear phase. 

Figures [7J and [8] show why in the high-temperature ferroelectric phase the polarisation 
vector is oriented along the b axis. Figure |7J gives a simplified view of the Mn layer in the ab 
plane, in which the spins inside squares depict the spins of Mn 4+ ions located inside oxygen 
octahedra, while the spins inside triangles are the spins of Mn 3+ ions in oxygen pyramids. 
The nearly collinear magnetic ordering in the high-temperature ferroelectric phase consists of 
antiferromagnetic chains along the a direction. 

The mechanism responsible for electric polarisation in this magnetic state involves, 
however, the ft! and ||| s P m chains along the b axis, such as shown in figure [8l These 
chains contain the polar Mn 4+ - Mn 3+ bonds, connecting parallel spins, and the Mn 3+ - 
Mn 4+ bonds with opposite polarity, connecting antiparallel spins. Importantly, the charge and 
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Figure 8. Electric polarisation induced by magnetostriction along the ^-chains. The two 
order parameters T7 1 and T72 describe degenerate ferroelectric states with opposite directions of 
electric polarisation. 

spin modulations in the chains have the same period, in which case the conventional exchange 
striction destroys the cancellation of electric dipoles of the polar bonds and induces electric 
polarisation along the chains, as illustrated in figure [8] 

This mechanism also works in the low-temperature (incommensurate) ferroelectric 
phase, which has the same periodicity in the ^-direction. The amplitude of the exchange 
striction is, however, the largest for collinear spins [see ©], which explains the drop of the 
polarisation at the transition to the low-temperature phase. For example, if the value of the 
magnetoelectric coupling oc f — j' 4 — J' 5 , is chosen such that the electric polarisation induced 
by the 'high-temperature' configuration shown in figure |6£ a) is 1000 /iC m~ 2 , then for the 
'low-temperature' configuration shown in figure Ob) it equals 500 /iC m~ 2 . 

4.3. Phenomenological approach 

In this subsection we discuss the phenomenological description of spin states and the 
magnetoelectric coupling mechanism discussed above, which will clarify the similarities 
between RM^Os an d other multiferroic materials. The positions of Mn 3+ and Mn 4+ in the 
paramagnetic unit cell are shown in tableland the eight vector order parameters can be found 
in table [2] (we use the notations of Bertaut et al 11351 ). 

Table 1. The coordinates of Mn ions, where x w 0.41, y w 0.35 and z ~ 0.26 (for BiMn 2 O s ). 
Mn 3+ Mn 4 + 



T\ = 


1/2) 


^5 


= (l/2,0,z) 


r 2 = 


(-*,-?, 1/2) 




= (1/2,0, -z) 


r3 = 


(1/2 -x, 1/2 +y, 1/2) 


m 


= (0,l/2,z) 


r 4 = 


(l/2+x,l/2-y,l/2) 


r 8 


= (0,l/2,-z) 



For discussion of phenomenological description of the magnetoelectric coupling, we 
have chosen the relatively simple case of BilVh^Os, which shows the commensurate spin 
ordering with Q = (1/2,0, 1/2). In this case the components of the order parameters belong 
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Table 2. Magnetic order parameters. 





Mn 3 + 




Mn 4 + 


F = 


S1 + S2 + S3 + S4 


F' 


= 5 5 + 5 6 + 5 7 + 5 8 


C = 


Si + 52 — S3 — S4 


a 


= 55 + 56 — 57 — 58 


G = 


Si — S2 + S3 — S4 


G' 


= 55 — 56 + 57 — 58 


A = 


Si — 52 — S3 + S4 


A 1 


= 55 — 56 — 57 + 5 8 



r 



r 2 



Table 3. Irreducible representations of the space group Pbam for Q = (1/2,0, 1 /2). 
2 X 2 y 2 Z m x m v m z I 

l\ /l 0\ /0 -A (\ 0\ (0 l\ (\ 0\ ( -1 

-A (\ o\ /o i\ /-1 o\ /o i\ (-1 o\ /o -1 

1 J V -1 J 1 J I 1 J v-i J l -1/ l-i 



Table 4. Basis vectors of the space group Pbam for Q = (1/2,0, 1/2). 




to one of the two two-dimensional representations, A or T2, of the Pbam group ll35l l48l . 
BiMn205 only shows the 'collinear' state with the a-components of the Mn 3+ and Mn 4+ 
spins described, respectively, by the order parameters F x = — 3.1/!b and G' x = 2.4/ie with 
small b components C y = — 0.8/ie and A' y = 0.6/ie ll48l corresponding to a small rotation 
between spins in neighbouring antiferromagnetic chains. 

f F x \ 

Since F x is the part two-dimensional representation I ^ I G A, the state described 

by the order parameter C x is another ground state of the system. These two states are 
related by inversion, which transforms F x into —C x and vice versa. It is easy to check that 
the magnetoelectric coupling of the form — X x P y (F^ — C^) is invariant upon all symmetry 
transformations of the paramagnetic phase, so that the order parameters F x and C x describe two 
ferroelectric states with opposite directions of electric polarisation. It also easy to check that 
the couplings —XyP y (Fy — C y ) and — X Z P Z {F} — Cf ) are also allowed by symmetry, which is 
a strong indication that the mechanism inducing electric polarisation in magnetically ordered 
state is invariant upon the global spin rotation and the coupling can be written in the form 
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Due to the exchange coupling between the Mn 3+ and Mn 4+ ions, the order parameters G' x 

A' 



between the two spin subsystems is phenomenologically described by — g (F X G' X — C X A' X ) ll20ll . 
Thus, more generally, the magnetoelectric coupling should be written in the form 

<t> me = -XPy(rf-r]i), (5) 

where f ^ 1 ) belongs to Ti representation. This form of the third-order magnetoelectric 

coupling was discussed previously in the context of the orthorombic manganites with the 
is-type magnetic ordering E2l and is typical for improper ferroelectrics ll25l . 

4.4. Static and dynamic dielectric susceptibility 

The contribution of the coupled spin-lattice degrees of freedom to static dielectric 
susceptibility is given by 



X\ 1 ~ 



hi U 



where 



Aij = J iJ (S i Sj)+8 ij 



£ K ia (l {Sihaf - Sf) - ^Jik (SiS k 



ce=l,2 



(7) 



and 



It^ftASixSj]^ (8) 

j 

The second term in © is the spin contribution to the dielectric constant due to virtual 
excitations of magnons by electric field (this will become more apparent in the discussion 
of the dynamic susceptibility). 

The last term in © describes the shift of the phonon frequency due to a change of the 
spring constants in magnetically ordered states. This effect is known in condensed matter 
spectroscopy as spin-phonon coupling Il4~9ll50ll . Phenomenologically, this effect is described 
by the fourth-order magnetoelectric coupling of the type P 2 L 2 , where L is a magnetic order 
parameter. In most cases, magnon and phonon branches coupled through this term experience 
'repulsion' and phonon hardens. However, in the magnetically frustrated compounds, this 
contribution to the dielectric constant can have either sign and can result in either hardening 
or softening of phonons in the magnetic phase (see, e. g., lf5TI on CdCr2S4 spinel). 

Equations of motion describing the coupled spin-lattice dynamics have the form, 

p _ Xo®Z dH 
r b — V dP h -> 

(9) 



S, 



H X S * 
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where Xo = X\^ + Xi an d 1S me t> are frequency of the polar phonon. Omitting the 
fourth order coupling term and solving linearized equations of motion, we obtain the dynamic 
dielectric susceptibility: 

hi 

J 0/ ' ij 



X- 1 (co) - Xo' 1 ( 1 - - ±J> \{BA- co 2 ) l B > 



where 



B U= J ij + % E (K ia (Sik a ) 2 -K i3 Sf)-^Y, J ik( S i S k)- (ID 
*i a=l,2 *i £ 

The second term in (fTOb describes the transfer of a part of electric dipole spectral 
weight from phonon to magnon frequencies, which turns magnons coupled to phonons into 
electromagnons. If such an electromagnon has lower frequency than the phonon, the static 
dielectric constant e(0) = 1 + 4nx(0) increases as a result of the coupling showing a step- 
like anomaly. Furthermore, frequencies of the mixed spin-lattice excitations (poles of the 
dielectric susceptibility (flOl) ) are shifted down with respect to the 'bare' magnon frequencies, 
found from 

det (BA - co 2 ) = 0. (12) 

Note that the electromagnon term in (flOl) disappears for collinear spin states, as /,-, defined by 
©, is zero in this case. This can be understood as follows. Classically, magnons correspond 
to spin oscillations that are orthogonal to ordered spin vectors. The change of the scalar 
product of a pair of collinear spins is then proportional to the second power of the amplitude 
of the oscillations. Since the magnetoelectric coupling in our model originates solely from 
Heisenberg exchange and only involves scalar products of spins, the linear coupling of electric 
field to magnons is absent for collinear spins and the lowest-order process is the photo- 
excitation of a pair of magnons. 

In figures |9t a) and (b), we plot the real and imaginary parts of the dielectric function (red 
and blue lines, respectively) for the 'high-temperature' and 'low-temperature' states shown 
in, respectively, figures [6^ a) and (b) (the imaginary part was obtained by the shift ft) — > ft) + /| 
with 7 = 1 K). As was explained above, the coupling of magnons to the electric component 
of light and significant electric dipole absorption at magnon frequencies is only present in 
the non-collinear 'low-temperature' phase, in agreement with experiment. This result may 
seem somewhat counterintuitive: while the spontaneous electric polarisation induced by the 
Heisenberg exchange striction is maximum for the collinear state, the excitation of magnons 
by the electric component of light (electromagnons) requires non-collinear spins and is only 
observed below the spin re-orientation transition. 

The red points in figure [9] mark the bare frequencies of the softest magnons with zero 
wave vector, found from (fT2~l) . Since the magnetic unit cell in this calculation contains 
16 mangetic ions, the total number of such magnons is also 16. However, only one of 
them is strongly coupled to the electric component of light and significantly contributes to 
the dielectric constant. This electromagnon corresponds to relative rotation of spins of the 
neighbouring antiferromagnetic chains, which gives rise to oscillations of the induced electric 
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Figure 9. (Colour online) The model calculation results: Frequency dependence of the real 
(El) and imaginary (£2) parts of the dielectric function (red and blue lines, respectively) for 
the 'high-temperature' collinear state (panel a) and the iow-temperature' non-collinear state 
(panel b) shown, respectively, in panels (a) and (b) of figure[6] Red points are selected magnon 
frequencies of spins decoupled from the lattice found as the roots of 1121 . A magnon at 
11.8 cm" 1 couples to a polar phonon at 100 cm -1 and becomes an electromagnon observable 
as the peak of the £2. 



polarisation in the b direction. The frequency of this uncoupled magnon for the non-collinear 
spin configuration shown in figure Ob) is 1 1.8 cm -1 . As the position of the absorption peak 
in figure Ob) is clearly lower that this magnon frequency, our parameters correspond to the 
strong magnetoelectric coupling case. The strong coupling is apparently necessary, if the large 
difference between dielectric constants of the 'high'- and 'low'-temperature phases Ae'(0) 
(= 3.25 in our calculation) is associated solely with the absorption peak emerging in the non- 
collinear state. 

The results of our model calculations presented in figure MX>) show that the symmetry 
properties of one magnon and the lattice allow an electromagnon in this system. Next, we 
discuss the strength of coupling. In our model, the first derivatives of the exchange integrals 7, 7 
are the coupling constants. Both magnetically induced polarisation and electromagnon were 
calculated using the same set of J' values. Characteristic value of /' used in our calculation 
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was dJ^/dy = 4 meV/A which is much less than, for example, 40 meV/A for ZnG^C^ [|52| . 
5. Discussion 

Our model calculations show that the observed values of the polarisation and electromagnon 
strength are obtained with realistic symmetric exchange constants. Also, the electromagnon 
coupling has been shown to be of Heisenberg type by the observed selection rules. These 
are strong arguments in favour of the Heisenberg exchange origin of the polarisation and 
electromagnons in 7?Mn205 family. However, in comparing the experimental and theoretical 
results, we should keep in mind that in the compounds of interest both symmetric and 
antisymmetric exchange mechanisms are operative, and they both may produce polarisation 
and electromagnons. Therefore, it is important to compare both the spontaneous polarisation 
and the electromagnon selection rules in both i?Mn03 and 7?Mn205 with these two 
mechanisms. 

In particular, the origin of ferroelectricity in i?Mn205 remains controversial. An 
alternative to the magnetostriction scenario discussed here is electric polarisation induced by 
the be spiral, as observed in recent experiments dl4"6ll4"7ll53l0 via the inverse Dzyaloshinskii- 
Moriya mechanism of relativistic origin ifTTl [T8l l24l . One natural question is whether the 
electromagnon study can clarify this controversy. 

Electromagnon excitations for the cycloidal spiral state were studied by Katsura et al 
|fl4"1|. According to their theory, the electric field of light that excites electromagnon has to 
be orthogonal to the direction of the electric polarisation induced by the spiral, which can be 
easily understood as follows. The spin spiral induces electric polarisation that lies in the spiral 
plane and is orthogonal to the propagation wave vector of the spiral ||24"1|. Thus the £>c-spiral 
in 7?Mn205 with the wave vector along the c axis would induce polarisation along the b axis, 
in agreement with experiment. An electric field applied in the direction perpendicular to the 
spiral plane will result in a small rotation of this plane, as a result of which the spontaneous 
polarisation vector acquires a component parallel to the applied field. An oscillating electric 
field orthogonal to the spiral plane will then excite oscillations of this plane, which is precisely 
the electromagnon of Katsura et al lfl4l . Since the polarisation lies in the spiral plane, the 
polarisation of light should be orthogonal to the spontaneous polarisation. Thus, to excite the 
oscillations of the be spiral would require e\\a, whereas in experiment electric field is parallel 
to the b axis and the direction of electric polarisation. Indeed, this can be an argument against 
the spiral scenario of multiferroicity of 7?Mn205. 

However, the situation is complicated by the fact that even in RMnOj, compounds, where 
the spiral origin of the magnetically-induced electric polarisation is well established, the 
selection rule e _L P is not obeyed. This was demonstrated by optical measurements in 
magnetic field. Magnetic field applied to DyMn03 OTI and TbMn03 [|43l gives rise to the 
magnetic flop transition at which the Z?c-spiral is replaced by the a^-spiral. The spiral flop, 
however, has no effect on the polarisation of electromagnons, which are always observed 
for e\\a. Also in Euo.75Yo.2sMn03, where the spin spiral lies in the ab plane already for 
zero magnetic field and the induced electric polarisation is parallel to a, the selection rule for 
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electromagnons remains the same: e\\a l|29ll . Therefore, the selection rules in both 7?Mn03 
and 7?Mn205 are tied to the lattice, not the spin plane. The selection rule for 7?Mn03 appears to 
originate from the GdFe03 distortions of the perovskite structure of orthorombic manganites, 
which generates new types of magnetoelectric interactions and couples light to magnetic 
excitations that are different from the oscillations of the spiral plane If54l . Furthermore, these 
distortions may also couple to the zone edge magnons and account for the broad peak observed 
in 7?Mn03 at high frequencies. However, an alternative mechanism for the broad peak is the 
photoexcitation of bi-magnons and so the precise origin of this high frequency peak in7?Mn03 
remains an open question. 

In any case, it appears that electromagnons in 7?Mn03 are also induced by 
symmetric exchange even though the spontaneous polarisation is produced by the DM-type 
antisymmetric exchange mechanism. This is understood as a consequence of the vanishing of 
any static polarisation induced by symmetric exchange due to the symmetry of the lattice lfT8l . 
Another important question then is why the antisymmetric exchange does not produce 
electromagnon excitations? We suspect that the relativistic interactions are just too weak 
to produce observable signals. The observation of these weak DM induced electromagnons 
would help clarify the overall picture of of spin-lattice interactions in multiferroics. 

Consequently, while our experimental and theoretical results strongly suggest that the 
magnetoelectric coupling in i^M^Os is governed by Heisenberg exchange, we cannot safely 
rule out the spiral scenario for the spontaneous polarisation in these materials on the basis 
of optical data. However, the results of recent neutron scattering experiments Il53ll suggest 
rather that the spiral components appear as a consequence of the existence of spontaneous 
polarisation. 

6. Summary 

We have presented experimental evidence and theoretical analysis that demonstrate that 
electromagnon excitations are present in both i^M^Os and 7?MnC>3 multiferroics. These 
electric dipole mixed magnon-phonon modes are observed in the far infrared and match 
excitations observed in inelastic neutron spectra in both classes of materials. The polarisation 
selection rules observed in the infrared experiments show that the electromagnon excitations 
are generated only by symmetric exchange in both classes of materials. The observed 
selection rules are tied to the lattice, not the spin plane, in contrast to the predictions of the 
antisymmetric exchange model. 

To theoretically account for electromagnons, we have considered third and fourth 
order coupling between the lattice and spins. The fourth order terms produce spin-phonon 
interactions that lead to shifts of the magnon and phonon frequencies near magnetic phase 
transitions observed in many magnetic materials. The third order coupling terms can produce 
mixed magnon-phonon excitations - the electromagnons - only for non-collinear spin orders. 
The third order terms are also responsible for the spontaneous polarisation in multiferroics J6]|. 
We also developed a simple microscopic model of 7?Mn205 that explains the transition 
between collinear and non-collinear spin states, the magnetically-induced electric polarisation 
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in the collinear state and the appearance of the electromagnon peak in the non-collinear 
state as well the polarisation of light for which this peak is observed. The mechanism of 
magnetoelectric coupling in this model relies entirely on isotropic Heisenberg exchange and 
magneto strictive coupling of spins to a polar phonon mode. 
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